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f"-) , Abstract. A class of special solutions are constructed in an intuitive way for the ultradis- 
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depending on the function-type and the system parameter. 
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1 Introduction 

Ultradiscretization pQ is a limiting procedure transforming a given difference equation into 
a cellular automaton, in which dependent variables also take discrete values. To apply this 
procedure, we first replace a dependent variable x n in the equation by 
> 

^2 : x n = e*»/«, (1) 

where e is a positive parameter. Next, we apply elog to both sides of the equation and take the 
t — ■ limit e — > +0. Then, using identity 

o. 

lim elog (e x/£ + e Y/e ) = max(X,Y), 

the original difference equation is approximated by a piecewise linear equation which can be 
regarded as a time evolution rule for a cellular automaton. In many examples, cellular automata 
obtained by this systematic method preserve the essential properties of the original equations, 
such as the qualitative behavior of exact solutions. However, the ansatz (Qj) is only possible if 
the variable x n is positive definite. This restriction is called 'negative problem'. 

From theoretical and application points of view, it is an interesting problem to study ultradis- 
crete analogs of special functions and their defining equations, including the Painleve equations. 
Ultradiscrete analogs for some of the Painleve equations and their special solutions are discussed, 
for example, in [21 El S] - However, the class of solutions for ultradiscrete Painleve equations has 
been restricted because of the negative problem. Some attempts resolving this problem are 
reported, for example, in [5j [61 [7]. The authors and coworkers study in [5] an ultradiscrete 
Painleve II equation with sinh ansatz and discuss its special solution of Bi function type. 

In order to overcome the negative problem, a new method 'ultradiscretization with parity 
variables' (p-ultradiscretization) is proposed in [Sj. The procedure keeps track of the sign of 
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original variables. By using this method, the authors and coworkers present [9] a p-ultradiscrete 
analog of the g-Painleve II equation (q-PII), 

(z(qr)z(r) + 1) {z(r)z{q~ l r) + l) = ^g-. (2) 

In [9] , we also discuss a series of special solutions corresponding to that of g-PII written in the 
determinants of size N. However, the resulting solutions are reduced to only one solution for the 
p-ultradiscrete Painleve II (udPII) equation. In this paper, we construct other series of special 
solutions for udPII and discuss their structure. In Section [2J we introduce the results in [9] for 
the p-ultradiscrete Airy equation. Then, we construct special solutions for udPII in Section [3l 
These solutions are, from their construction, considered to be counterparts of those of g-PII 
written by the determinants. Finally, concluding remarks are given in Section [U 

2 Ultradiscrete Airy equation with parity variables 

We start with a g-difference analog of the Airy equation 

w(qr) — tw(t) + w[q~ r) = 0, (3) 

which reduces to the Airy equation 

d 2 v 



dsi +SV = ° 

in a continuous limit. 

In order to ultradiscretize (J3]), we put r = q m and q = e®' £ (Q < 0). Furthermore, we 
introduce an ansatz for p-ultradiscretization, 

w(q m ) = { 8 (u: m )-8(- Um )}e Wm/e , 

where uj m € {+1, —1} denotes the sign of w{q m ) and s(u) is defined by 

Taking the ultradicrete limit, we obtain a p- ultradiscrete analog of the Airy equation 

max (W m+ i + S(uj m+1 ),mQ + W m + S(-u m ), W m -\ + S(u ro _i)) 

= max(M / m+ i + S(-uj m+1 ),mQ + W m + S(u) m ),W m -i + 5(-w m _ 1 )), (4) 

where S(u) is defined by 
S(u) 





An ultradiscretized variable is represented by a pair of u m and W m , which is denoted as W n = 
(uj m ,W m ) in what follows. It is possible to rewrite the implicit form @ into explicit forward 
schemes 

7. 1 n sgn(F m ), co m = -w m _i or F m ± 0, 

indefinite, oj m = w m _i and F m = 0, 
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W, 



m+1 



max (mQ + W m , W m -i) 



<W, 



m— 1) 



uj m = -Wm-i or F m + 1 0, 
L0 m = w m _i and F m = 0, 



where F m := mQ + W m — W m —\. Note that we generally have both of unique and indeterminate 
schemes depending on given values of (io m ,W m ) and (w m _i, W m _i). The explicit backward 
schemes are obtained by replacing m ± 1 with m+" 1, respectively. 

We find two typical solutions of (J1J). One is an Ai-function-type solution for Wo = (+1,0) 
and Wi = (+1,0), 



uAi(m) = (u} m ,W m ) 



( ( , , m(m-l) 

(-i) 2 .0,, 



+1 m(m-l) Q 



m > 0, 
m < — 1, 



and the other is a Bi-function-type for Wo = (+1,0) and Wi = (—1,0), 

m > 0, 

77i(?n + 1) 



uBi(m) = (u; m ,W n 



m ( 77i + 1 ) \ 

(-i)-4-\ol, 



+i, 



Q , m<-l. 



They show similar behavior as those of the Ai and Bi functions, respectively. 

3 Ultradiscrete Painleve II equation with parity variables 

For the following discussion, we first introduce the results for (J2J) . It has been shown in |10| that 

r s (A %)<? (iv+i v) 



*W(r) 



q N g( N )(qT)g( N +V(T) 

g (N) {T)g{ N + l) {qT) 



N>0, 
N <0 



(5) 



^JV+l 9 (7V)^ T ) 5 (7V+l)( r )' 

solves ([2]) with a = q 2N+1 , where the functions g( N >(t) (N € Z) satisfy the bilinear equations 



,2 V A,+1) (g -l T)ff (AT) (g 2 r) _ q N Tg (N + l) {T)g (N } {qr) + gi N + l) ^(N) (r) = Q > 

W+i-Hq-^gWfr) - q 2N rg( N ^(r)g( N Hr)+ g ( N ^\qr) g ( N H q ^r) = 

for iV > and 

g 2JV + 2 5 ( JV+1 )(g- 1 r) 5 W(g 2 r) - g^+W"* 1 ) (r)</W (gr) + <? (7V+1 V)^(r) = 0, 
^ +2 5 (JV+1) (g- 1 r) 5 W(gr) - q 2N+1 rg( N+l Hr)9 {N) (r) + g^ +l \qr)g^{q-\) = 



(6) 

(7) 

(8) 

0) 



for iV < 0. It is also known that g^ N '(T) are written in terms of the Casorati determinant of 
size \N\ whose elements are represented by the solutions of ([3]). 

In order to construct ultradiscrete analogs of these equations, we put r = q m , q = e®' £ (Q < 0) 
and a = e ' £ . Furthermore, we introduce 

Then ([2]) is reduced to udPII, 

max Z m+1 + 3Z m + Z m _i + max{5(C m +i) + S(( m ) + S(C m -i), 
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S(-( m+1 ) + S(Cm) + 5(-Cm-l), S(-Cm+l) + S(-( m ) + S(Cm-l), 
S((m+l) + S(— Cm) + S(— Cm-l)}) -^m+1 + 2Z m + 5(Cm+l), 
2Z m + Z m _i + 5 (Cm-l), Z m + S(C m ), Z m + A + 2mQ + S'(Cm), 

Z m+ i + 2Z m + Z m _i + mQ + max{5'(-Cm+i) + <S(Cm-i), #(Cm+i) + ^(-Cm-l)}, 
Z m+ i + Z m + mQ + max{5(-Cm+i) + 'S'(Cm), S(Cm+l) + S(-Cm)}, 
Z m + Z m _i + m(5 + max{5(-Cm) + S'(Cm-l), ^(Cm) +<S(-Cm-l)} 
Z m+ i + 3Z m + Z m _i + max{5'(-Cm+i) + S'(-Cm) + S(-Cm-i), 

5(Cm+l) + 5(-Cm) + 5(Cm-l), 5(Cm+l) + S(Cm) + S(-Cm-l), 
^X - Cm+l) + S'(Cm) + £(Cm-l)})^m+l + 2Z m + S( — Cm+l] 



For 



and 



max 



max 



■ T J-/ ' — \jl IV / • — \15 1 1 o — J. / I j — (/t^_i_ ■ lib • — \_ -3 1 1 (,t J- / " 

2Z m + Z m _i + 5(-Cm-l), Z m + 5(-Cm), 2,71 + A + 2mQ + 5(-Cm), 

Z m+ i + 2Z m + Z m _i +mQ + max{S'(C m +i) + S(Cm-i), S(-( m+1 ) + S(-Cm-i)}, 

Z m+ i + Z m + mQ + max{5(C m +i) + S(( m ), S(-Cm+i) + S(-(m)}, 

Z m + Z m _i + mQ + max{5(Cm) + S(Cm-i), ^(~Cm) + 5(-Cm-i)}, "i,<2 • (10) 

and ((7J), we have their ultradiscrete analogs 

}NQ + G^ + G™ + max {s( 7 ™) + S( 7 $ a ) , SH^) + 5(- 7 ^ 2 ) } , 

+ max {5( 7 ^ +1) ) + S(-J ( S) • ^(-7^ +1) ) + SfrSi) } , 

GS 1} + Gg° + max {5(7^ } ) + S( 7 ^) ,S(-^ + + i ) )+S(- 1 ^)\ 

[2NQ+ G^ + G™ + -X {5( 7 S 1} ) + S(-^) , S(-^) + 5( 7 ffl 2 ) } , 
(A^ + HQ + G^+G^ 

+ max{5(7m V+1) ) +S(^ 1 ),S(-r^ +1) ) + SH*2i)}. 

^ +^> +max{s(^ (11) 



= max 



max 



max 



2NQ + GS 1} + G™ + max {s^) + S(^) , S(- 7m !^) + S^ift) } 

(2A r + m)Q + G^ +1) + G^> 

+ max{5(7^ +1) )+5(- 7 W),5(- 7 ( ) f +1 )) + 5( 7 ( ) f))}, 

G2£ 1} + Gfi2x + max {s( 7 £j 1} ) + 5( 7 Sx) - S(-^) + S(- 7 S) } 

+ max {5( 7 iS 1} ) + S(-7£2i) ■ SHL^) + stfi) } , 

(2N + ro )Q + GL W+1) + G™ 

+ max{5( 7 ^ +1 ))+5( 7 W),5(- 7 ^ +1 ))+S(- 7 W)}, 

^ + Gfi2i + max {s^ 1 >) + ^a.^ffiVtfi)} 



(12) 
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respectively. For and ©, we have 



max 



2(N + 1)Q + G^ + G iN) 



m+2 



■N + m + VQ + G^+G^, 

S(^ +1) ) + 5(- 7 if + \),5(- 7 Cf+ 1 )) + 5(7^0} , 



+ max 



G^ + QW + max {5( 7 SJ 1} ) + ^7^) > ^ffi 15 ) + *H£°) } 



max 



2(n + i)q + g^ 4 



(N) 

m+2 



+ max {S( 7 W>) + 5(- 7 S 2 ) , 5(- 7 S5 1) ) + %S 2 ) } , 

(iV + m + l^ + G^+G^ 

+ max {5( 7 ^) +5( 7 iSi). 5(-7^ +1) ) +5(- 7 ffi 1 )}, 

G2S 1} + GS° + -ax {5( 7 £S 1} ) + 5(-7^) , ShffiVW 



and 



max 



2(^+i)o + g2L + i 1) + g£J 1 



+ max{5( 7 S 1 )) +%2? 1 )^HS5 1) ) +^(-7^)}, 

(2AT + m + 1)Q + G^ +1 ) + GW 

+ max {S( 7 ^ +1) ) + S(- 7 W), S(- 7 ^ +1) ) + S( 7 M) } , 



g!S 1} + cffi m 



max 
+ max 



2(iV + 1)Q + G^jp -4 



(N) 

m+1 



(2N + m + 1)Q + G^ +1 ) + gL w > 

+ max{5( 7 ( / f +1 ))+5(7^),5(-7(f +1 ))+5(-7^ ) )}' 



G 



(JV+1) 
m+l 



+ G 



(TV) 
m— 1 



+ max 



respectively. Finally, the transformations ([5]) are reduced to 



AN) 

Sm 



(AT) (JV+l) (JV+1) (N) 
Im I m+l Im /m+l' 



ZL N) = GW + G <™> - G^ ^ - G™ - iVQ 



for iV > and 



f(«) = ^(N)JN+l) (N+l) (N) 
Sm /m /m+l Ira /m+l' 

4f ) = GW + G™ - g^ +1 > - cffi - (iV + 1) 



(13) 



H7sv^hf\)^hr)+s(a}], a*) 



(15) 
(16) 

(17) 
(18) 



for iV < 0. If we find solutions for the ultradiscrete bilinear equations, special solutions for 
udPII are obtained through (fT5j) - (fT8|) . 
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Hereafter we consider only the case of A = (2N + 1)Q in (|10p . which corresponds to a = q 2N+1 
in the discrete system. Firstly, we present the results reported in [5], that is, the Ai- function- 



type solutions for N > 0. Solutions of (fl"3j) and (fl~4"j) are given by Q m y ' = ("fm'fGrn') for 



,(JV) 



(N) n (N), 



N = 0,1,2, 



Im 



where 



' (N), i ., ™(m-l)N 

7o (-1) 2 



m > 0, 
m < — 1, 



gW 



' mJV(Ar - 1) 
2 

mN(m + iV 



q + g 



(AT) 
' 



-Q + G 



(TV) 
' 



m > 0, 
m < -1. 



Since Gm — >• — oo as m — > — oo in the same way as the uAi function, we call these solutions 
the Ai-function-type solutions. From these solutions, we have only one special solution of udPII 
with A = (2iV + 1)Q for N = 0, 1, 2, . . . 



Z (N) (AN) Z (N)\ 



((-1)"\0), m>0, 

( + 1,771(5) ) m < ~~ 1) 



(19) 



which does not depend on N. We note that Z 



(N) 



oo as m 



-oo. 



Secondly, we investigate Bi-function-type solutions for N > 0. We find that Q 



(o) 



and £. 



(i) 



(+1,0) 



uBi(m) solve (jlip and (|12p with AT = 0. By using this result, we inductively 



construct solutions Q. 



(AT+l) 



of the equations with Af > 1 for a given function Q m ' and assigned 



values of Qq and Q\ . We further assume that Q\ and Q\ are chosen so that 

Qm for any m are uniquely determined in (jlip and (|12p . Then we have the following solutions 
Gin = (7m , G m ), where 



7, 



(TV) 



(-1) 2 7o 7 

(m-2)(m-l)m(m + l) (jy) 

(-1) 24 70 , 

m(m+l)(m + 2)(m+3) ( M\ 

(-1) 24 7o , 

iv(iv-i) fjV) 
2 7o ; , 



(-1) 

rmN(N 



m > —1, 

-2 > m > -2N - 1, N : even, 
-2 > m > -2N - 1, N : odd, 
m < -2N - 2, 



G (JV) 



2 

mN(N 



^Q + Gi N \ 



2 

mN(N 



i ) Q+ (m-2M2m + l) Q + g ( JV); 



24 



1 



Q 



771 



l)(m + l)(2m-3) G (jv) 



24 
mAT(m + AT) Q _ JV(Ar-l)(4JV + l) Q + ^ 



m > —1, 

-2,-4,.. 
-3,-5,.. 

m < -2 A" - 2. 



/// 



/» 



c- /o(AT) 

hmce Gm. 



-2iV, 
-2AT-1, 



2 6 

-> — oo in the same way as the uBi function, we call these solutions the 



oo as m — > - 
Bi-function-type solutions. 

By substituting these solutions into (|15p and f)16[) . we obtain special solutions of udPII, 



Z (JV) = (AN) Z (N)\ 

m ISm > m J 



'((_l)m-i 9 0), m>-2N-l, 

(+l,-(m + 2AT + l)Q), m<-2AT-2. 



(20) 
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We notice that (|19p and (|20p have different asymptotic behavior in Z m for m — > — oo and in 

(0) 



the phases for m > 0. Furthermore, we remark that (|20p have iV-dependence. 
Thirdly, we study Ai-function-type solutions for N < 0. We find that Q. 



>(-!) 



(+1,0) and 
uAi(?n — 1) solve (j!3j) and (J14D with iV = — 1. Starting from these simple solutions, we 



inductively find the solutions Q m = (j m ,G m ) for N = 0, — 1, —2, . . . , where 



im 



G^ 



■ ( m+J v)( m+J v_i)iv , N s 
(-1) 2 <y\"\ m > -AT, 



(AT) 
.71 ' 



m < -JV, 



{ N(N + l)(m + N-l) „ {N) 
(^ + Cr x , 

-JV(m + iV-l)(m-l) ^ (TV) 
U> + Li 1 , 



m > -iV, 
m < -N. 



Substituting these Qm into (|17p and (|18p . we have special solutions of udPII with A 
(2N + \)Q for N 



-1,-2,... 



Z (iV) = (AN) Z (N)\ 



((_l)m-l >0 ), m>-N, 

(+l,-(m + 2N + l)Q), m<-N-l. 



(21) 



Typical behavior of these solutions is shown in Fig. [TJ They converge to as m — > — oo and 
oscillate for m > —N. It is interesting to note that (|21|) constructed from the uAi function has 
essentially the same structure as (|20p constructed from the uBi function. 



f;r"ex P (z;r") 





Figure 1. Behavior of special solutions (pH]) with Q 
respectively. 



-1. (a) and (6) are for TV = -1 and N 



Finally, we study Bi-function-type solutions for N < 0. We find that Q m ' = (+1,0) and 
7m = uBi(m — 1) solve (fT3j) and flT3J) with N = — 1. We construct solutions Q m of the 



(A+l) 



,(N) 



,(N) 



We 



equations with A^ < — 2 for a given function Q m and assigned values of Q\ and C/ : 
further assume that Q[ and ^ ' are chosen so that Q m for m < are uniquely determined 
in (fT3l) and (fT4|l . We again inductively obtain the solutions Q m = (7m , G m ) , where 

m > 3 — N, N : even, 

m > 3 - N, N : odd, 
N -l<m<2- N, 
m < N -2, 
(m-N - 3)(m - iV - 4)(m -N- 5)(m - iV - 6) 






(-1) 8 7i > 

( m+i y-3)(m+iV-4) (iV-l)(iV-3) -. fjyN 
(-1) 2 + 8 +i 7 } ; , 

(_l)PW,m 7 W > 

, , . (JV+2)(iV+3)(iV+4)(iV+5) , (-JV-, 
(-1) 24 +^} ', 



PiV,: 



24 



+ 



(N + 2){N + 3)(N + A){N + 5) 



24 
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( (m-l)N(N + l) „ N(N + 2)(2N-1) „ {N) 

2 Q 24 g + Gl ' 

m > 3 — N, N : even, 
( m -l)iV(iV + l) (JV-l)(iV + l)(2iV + 3) {N) 

2 g 24 g+ x ' 

m>3-N,N : odd, 
(m - l)iV(m + 3iV + 2) (m - 2)m(2m + 1) „(/v) 
4 Q + 24 <y + W , 



g™ 



N — l<m<2 — N,m: even, TV : even, 

(m - l)N(m + 3iV + 2) (m - 2)m(2m + 1) + 6 (JV ) 

4 ^ + ^4 g + Cri , 

N— l<m<2 — N,m: even, iV : odd, 

(m - l)iV(m + 3iV + 2) (m - l)(m + l)(2m - 3) (JV) 

-t/ H — v + "i ! 



24 
N -1 <m<2- N,m: odd, 



miV(m + AH iV(2iV 2 - 15iV _ 14) (JV) 



24 



m < N — 2,N : even, 



mAT(m + iV) (JV + 1) (2iV 2 _ 17JV + 3) 



-Q + 



24 



Q + Gi 



(TV) 



m < N -2,N : odd. 
From these solutions, we have special solutions of udPII, 

Z (N) = (AN) Z (N)\ 



((-i) m ,o), 

m-N 171 + N 



(-iy 



Q 



m-N-i m + N + 1 
(-1) 2 , 2 Q 

(+l,mQ), 



m > -iV- 1, 

|m| < — iV — 2,m — N : even, 

[ m | < _j\r -2,m-N : odd, 

m < iV + 1. 



(22) 



Typical behavior of these solutions is shown in Fig. [2j Note that (|2ip and ([22]) have different 
asymptotic behavior in Zm as m — > —00 and in the phases for m > — iV. We also comment 



that, although (122p is similar to (|19p . (122p has more complicated internal structure. 



fir'Wizlr") 



i'^xptzir"') 





Figure 2. Behavior of special solutions 
respectively. 



(6) 
with Q = —1. (a) and (6) are for TV = — 1 and N = —6, 
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4 Concluding remarks 

In this paper we have presented a class of special solutions for the p-ultradiscrete analog of q-PII. 
The solutions are classified into four groups; Ai-function-type and Bi-function-type solutions for 
the system parameter N > 0, and those for N < 0. In the preceding paper [9] are given only 
the Ai-function-type solutions for N > 0, which do not depend on N. Three other groups which 
are newly given in this paper do depend on N. Moreover, the solutions of each group have 
different structures. For example, we observe differences between the Ai- and Bi-function-type 
solutions in their asymptotic amplitude and phases, which may reflect the structure of solutions 
of difference and continuous equations. The Bi-function-type solutions for N < have fairly 
complicated internal structure, although we do not know the origin of these structures yet. At 
any rate, these results may indicate the richness of solution space of the ultradiscrete equation. 

For the continuous and discrete Airy equations, linear combination of Ai and Bi functions give 
their general solutions. In the ultradiscrete case, max(/, g) corresponds to the linear combination 
of functions / and g. Hence, we believe that the cases we treated in this paper cover quite wide 
class of special solutions of the ultradiscrete equations. 

Our method of constructing solutions is intuitive and purely based on the ultradiscrete equa- 
tions. We believe that the solutions we obtain correspond to those of q-PIl represented by the 
Casorati determinant of size \N\ whose elements are given by the q-difference Ai or Bi function. 
It is a future problem to clarify the relationship between discrete and ultradiscrete solutions 
through a limiting procedure. It is also a future problem to construct p-ultradiscrete analogs of 
other Painleve equations and their special solutions. 
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